Optimization of Network Robustness to Waves of Targeted and Random Attacks 
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We study the robustness of complex networks to multiple waves of simultaneous (i) targeted 
attacks in which the highest degree nodes are removed and (ii) random attacks (or failures) in which 
fractions pt and pr respectively of the nodes are removed until the network collapses. We find that 
the network design which optimizes network robustness has a bimodal degree distribution, with a 
fraction r of the nodes having degree fe2 = ((fc) — 1 + r)/r and the remainder of the nodes having 
degree fci = 1, where (fc) is the average degree of all the nodes. We find that the optimal value of r 
is of the order of Pt/Pr for Pt/Pr <C 1. 
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PACS numbers: 89.20.Hh, 02.50.Cw, 64.60.Ak 

Recently, there has been much interest in the resilience 
of real-world networks to random attacks or to attacks 
targeted on the highest degree nodes [1I1I1I1IEI1I1I3. 
Many real-world networks are robust to random attacks 
but vulnerable to targeted attacks. It is important to 
understand how to design networks which are optimally 
robust against both types of attacks, with examples be- 
ing terrorist attacks on physical networks and attacks by 
hackers on computer networks. Studies to date |3,|a| have 
considered only the case in which there was only one type 
of attack on a given network — that is, the network was 
subject to either a random attack or to a targeted attack 
but not subject to different types of attack simultane- 
ously. 

A more realistic scenario is one in which a network is 
subjected to simultaneous targeted and random attacks. 
This scenario can be modeled as a sequence of "waves" 
of targeted and random attacks which remove fractions 
Pt and Pr of the original nodes, respectively. The ratio 
Pt/Pr is kept constant while the individual fractions pt 
and Pr approach zero. After some time (after m waves of 
random and targeted attacks) the network will become 
disconnected; at this point a fraction fc — m{pt + Pr) of 
the nodes have been removed. This fc characterizes the 
network robustness. The larger fc, the more robust the 
network is. We propose in this Letter a mathematical 
approach to study such simultaneous attacks and find 
the optimal network design one which maximizes fc In 
our optimization analysis, we compare the robustness of 
networks which have the same "cost" of construction and 
maintenance, where we define cost to be proportional to 
the average degree (fc) of all the nodes in the network. 

We study mainly two types of random networks: 

(i) Scale-free networks. Many real world computer, so- 
cial, biological and other types of networks have 
been found to be scale free, i.e., they exhibit de gree 
distributions of the form P(fc) ~ fc"^ H Ellil 
EIEIEIEEESIII- For large scale-free networks 



with exponent A less than 3, for random attacks es- 
sentially all nodes must be removed for the network 
to become disconnected jj,!^. On the other hand, 
because the scale free distribution has a long power- 
law tail (i.e. hubs with large degree), the scale- free 
networks are very vulnerable with respect to tar- 
geted attack. 

(ii) Networks with bimodal degree distributions. For re- 
silience to single random or single targeted attacks, 
certain bimodal distributions are superior to any 
other network |3, Ig • Here we ask if these networks 
are also most resilient to multiple waves of both 
random and targeted attacks. 

We present the following argument that suggests that 
the degree distribution which optimizes fc is a bimodal 
distribution in which a fraction r of the nodes has degree 



fc2 = 



(fc) -1-Hr 



(1) 



and the remainder has degree fci = 1 and we show that 
r is of the order oipt/pr. To optimize against random 
removal, we maximize the quantity k, = (fc^)/(fc), since 
for random removal the threshold is yi 
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Since we keep (fc) fixed, k, is just the variance of the de- 
gree distribution and is maximized for a bimodal distri- 
bution in which the lower degree nodes have the smallest 
possible degree fci = 1 and the higher degree nodes have 
the highest possible degree consistent with keeping (fc) 
fixed, fc2 = ((fc) — 1 -l- r)/r. Thus, fc2 is maximized when 
r assumes its smallest possible value, r = 1/N. On the 
other hand, if all of the high degree nodes are removed 
by targeted attacks, the network will be very vulnerable 
to random attack. So we want to delay as long as pos- 
sible the situation in which all of the high degree nodes 



are removed by targeted attacks — which argues for not 
choosing r as smah as possible but choosing r such that 
some high connectivity nodes remain as long as there 
are some low connectivity nodes. Such a condition is 
achieved when r is of the order oi pt/pr- 
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FIG. 1: (a) The threshold fc of three bimodal networks with 
(k) = 3, with (i) r = 2 X 10"^ and fe = 200, (ii) r = 5 x 10"^ 
and fc2 ~ 90, and (iii) r = 10"'^ and fe = 50. The results are 
plotted as a function of the ratio Pt/Pr for three fixed values 
of Pr. These plots show that the values of the threshold are 
dependent only on the ratio pt/Pr and independent of the 
value of Pr itself, (b) Scaled plot of the data in (a) . The data 
show that the plots collapse in the region where r/ (pt/pr) ^ 
1. 

The method we employ for determining the threshold 
makes use of the following: the general condition for a 
random network to be globally connected is 



(fc) 



> 2. 



(3) 



Random removal of a fraction pr of nodes from a network 
with degree distribution Po{k) results in a new degree 
distribution |3 
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where K is the upper cutoff of the degree distribution. 
Targeted removal of a fraction pt of the highest degree 



nodes reduces the value of upper cutoff K to K, which 
is implicitly determined by the equation 
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(5) 



k=K 



The removal of high degree nodes causes another effect. 
Since the links that lead to removed nodes are elimi- 
nated, the degree distribution also changes. This effect 
is equivalent to the random removal of a fraction of p 
nodes where 
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The average (fc)o is taken over the degree distribution 
before the removal of nodes J5| . Equation Q with pr re- 
placed by p can then be used to calculate the effect of the 
link removal. Starting with a certain initial degree distri- 
bution, we recursively calculate P{k) alternating between 
random and targeted attack using Eqs. I@J, IJSJ), and ||HJ), 
and calculate k after each wave of attacks. When k < 2 
global connectivity is lost and fc — m{pr + pt) where m 
is the number of waves of attacks performed. 

We begin our study by first establishing numerically 
that, for small values of pt, pr and pt/pr, the thresh- 
old fc depends only on pt/pr- In Fig. QJa), we plot the 
threshold fc of a network with a bimodal degree distri- 
bution with (fc) = 3 for various values of Pr and r as 
a function of the ratio pt/Pr- The collapse of the plots 
with the same r but different Pr shows that the values of 
threshold are essentially independent of the value of Pr 
itself but depend only on the ratio pt/Pr u^- 

In Fig. ^b) we plot fc against the scaled variable 
r/{pt/pr)- We see that the plots for different values 
of r collapse, indica ting that only the scaled variable 
r/{pt/pr) is relevant [2J|. 

Next we study the dependence of fc on fc2 • As seen in 
Fig. 121 as expected the maximum values of fc for various 
values oipt/pr are obtained when fc2 is maximum (i.e., 
when fci = 1, see Eq. Q.). 

We are now in a position to determine the value of r 
which optimizes fc, ^opt. In Fig.|3 we plot fc as a func- 
tion of the scaled parameter r/{ptlpr) with fc2 set to the 
maximum value possible for each value of r . We note that 
there is a transition at a well-defined value of r / {pt/pr) 
at which fc increases rapidly to a shallow maximum /°p* 
at ropt/{pt/Pr) « 1.7. This value oi ropt/iPt/Pr) is valid 
for pt/pr ^1. In order to determine ?'opt /(pt/Pr) over 
a wider range, we make extensive numerical calculations 
for 10"'^ < pt/Pr < 0.1. For each value oi pt/pr, we cal- 
culate the value ''opt /(pt/Pr) where fc takes its maximum 
value and find 
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FIG. 2: The threshold /c versus ^2 for a bimodal network with 
(fc) = 3 and r = lO"'^ for three values oipt/Pr- The value of pr 
is fixed at 0.02. For each value oi pt/pr, the thresholds take 
their maximum values at the maximum k^ (obtained when 
fci = 1). 
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FIG. 3: The threshold fc versus the scaled parameter 
r/{pt/pr) for a bimodal network with (fc) — 3 and ^2 maxi- 
mum (i.e., fci = 1). 



within the range of our calculation. For larger values of 
Pt/Pr, '"opt = 1 and from Eq. |^ all nodes have degree 
(k). In Fig.^ we plot the values of the optimal threshold 
jopt i^y ,^ thick solid curve. 

In Fig. 21 we also plot the values of the threshold fc for 
the same bimodal network but we fix r independent of 
Pt/Pr- We see that these configurations are not signifi- 
cantly less robust than the optimal configuration. Thus, 
even if we do not know the ratio pt/pr exactly we can 
design networks which will be relatively robust. For ex- 
ample, the bimodal network with r = 0.03 is relatively 
robust for pt/pr ^0.1 and the bimodal network with 
r = 0.09 is robust for pt/pr ^ 1- Also plotted in Fig. ^ 
is the optimal scale- free network with {k) = 3. We see 
that the optimal bimodal network is more robust than 
the optimal scale-free network and we can even pick a 
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FIG. 4: The threshold fc versus pt/pr- The topmost (thick- 
est) curve is for a bimodal network with (fc) = 3 with 
fci = 1 and with r optimized by Eq. ^ for each value 
of Pt/Pr- The values of the threshold for the same bi- 
modal network with fci = 1 when we fix r independent 
of pt/pr are plotted in thin curves. The values of r are 
r = 0.001,0.002,0.005,0.01,0.03,0.05,0.07,0.09,0.11,0.13, 
and 0.15, from left to right. The curve marked with crossed 
circles (©) is a plot of the threshold values for a scale-free net- 
work with (fc) = 3, A*' = 10'', and with exponent chosen for 
each pt/Pr to optimize the threshold. Note that the thresh- 
olds for bimodal networks with 0.03 < r < 0.09 are always 
more robust than the optimized scale-free network. 



configuration with fixed r (e.g. r = 0.03) which is more 
robust than the optimal scale-free network in most ranges 
oipt/pr- 




FIG. 5; Realization of bimodal network with A*' — 100 nodes, 
(fc) = 2.1 and r = 0.1, so there are rN — 10 "hub" nodes of 
degree 12, as found from Eq. 10. 

In Fig. |5l we show a typical optimal realization of a 
bimodal network. The network of iV = 100 nodes con- 
sists of rN nodes with k — k2 { "hubs" ) which are highly 
connected among themselves; the nodes of single degree 
are each connected to one of these hubs. We note that 



while the hubs are highly connected among themselves, 
they do not form a complete graph — every hub is not 
connected to every other hub. For larger A^, the fraction 
of hubs to which a given hub connects decreases but the 
robustness of the network is unchanged. 

In summary, we have provided a qualitative argument 
and numerical results which indicate that the most ro- 
bust network to multiple waves of targeted and random 
attacks has a bimodal degree distribution with a fraction 
r of the nodes having degree fe — ((^) ~ 1 + r)/r and 
the remainder of the nodes having degree 1. The opti- 
mal value of r is approximately 1.7 (pt/Pr) for Pt/Pr- ^ 1. 
For larger values oi pt/pr, the optimal value of r is 1 and 
all nodes have degree (fc). Even ii pt/pr is not known 
exactly, a value of r can be chosen which maximizes the 
network robustness over a wide range of values oi pt/pr, 
as seen in Fig. 0| 

We note that while the optimal distribution found here 
and that found in Ref. jSi] are both bimodal, the values of 
the parameters characterizing these distributions are dif- 
ferent. As found in Ref. [S| the network with optimal re- 
silience to either random or targeted attack has r = 1/N 
and k2 ~ r~'^l^ . Finally, we note that it is possible to 
prove analytically that for the case in which a single tar- 
geted attack followed by a single random attack results 
in the network becoming disconnected, the optimal dis- 
tribution is also bimodal with k\ = \,ki^ {{k) — l + r)/r 
and r of the order oi pt/pr [25J — supporting the results 
found here for multiple waves of attacks. 
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